Exact solution of a spin-1/2 XX chain with three-site interactions in a random 
transverse field: Influence of randomness on quantum phase transition 
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We present exact results for the ground-state and thermodynamic properties of the spin-1/2 XX 
chain with three-site interactions in a random (Lorentzian) transverse field. We discuss the infiuence 
of randomness on the quantum critical behavior known to be present in the nonrandom model. We 
find that at zero temperature the characteristic features of the quantum phase transition, such as 
kinks in the magnetization versus field curve, are smeared out by randomness. However, at low but 
finite temperatures signatures of the quantum critical behavior are preserved if the randomness is 
not too large. Even the quantum critical region may be slightly enlarged for very weak randomness. 
In addition to the exact results for Lorentzian randomness we present a more general discussion 
of an arbitrarily random transverse magnetic field based on the inspection of the moments of the 
density of states. 

PACS numbers: 75.10.Jm 

Keywords: quantum phase transitions, random quantum spin chains, multi-site interactions. Green functions, 
density of states 



I. INTRODUCTORY REMARKS 

In recent years the theory of quantum phase transi- 
tions has been in the focus of very active researchji"— 
The quantum phase transitions take place at zero tem- 
perature by changing a control parameter and emerge as 
a result of competing different ground-state phases. Im- 
portantly, quantum phase transitions can influence the 
behavior of systems over a wide range of the phase dia- 
gram at nonzero (sometimes quite large) temperatures. 
Exactly solvable quantum models exhibiting a quantum 
phase transition are notoriously rare. A well-known ex- 
ample of a solvable model is the spin-1/2 Ising chain in 
a transverse field, where a zero-temperature transition 
from the ordered quantum Ising phase (small transverse 
fields) to the disordered quantum paramagnetic phase 
(large transverse fields) takes place. This model is of- 
ten used for illustration of basic concepts in the quan- 
tum phase transition theoryi^"— In general, spin-1/2 XY 
chains^ provide an excellent ground for various statistical 
mechanics studies since in many cases the calculations 
can be performed without any approximation. More- 
over, there are some real-life compounds which can be 
viewed as realizations of one-dimensional spin-1/2 XY 

Quite recently, two other classes of solvable mod- 
els have been found, namely a two-dimensional Ki- 
taev modeli^ and a spin-1/2 XY chain with multi-site 
interactions^^"— (see also Ref. [isl ). The model belong- 
ing to the latter class is of interest in this paper. This 
model has an essentially richer ground-state phase di- 
agram as the standard one-dimensional spin-1/2 XY 
model. In particular, it may exhibit several gapless spin- 



liquid phases and quantum phase transitions between 
them^^ 

On the other hand, quantum models with random 
Hamiltonian parameters present another class of mod- 
els for which an exact solution cannot be found easily. 
A solvable model with diagonal Lorentzian disorder was 
introduced by LloydJ^ Later on Lloyd's idea was used 
to study random spin-1/2 XX chains»22, Also an exten- 
sion to correlated off-diagonal Lorentzian disorder and 
its application to spin-1/2 XX chains was considered, 
see Refs.[2l|[23. 

Naturally, the investigation of quantum phase transi- 
tions in systems with randomness is a challenging task. 
The random transverse-field Ising spin chain is known as 
a tractable model to study effects of quenched random- 
ness on critical behavior Within the context of random 
quantum systems exactly solvable models may play an 
important role. Merging together the above mentioned 
solvable quantum spin models with three-site interactions 
and Lloyd's model of disorder we present here an exact 
analysis of a specific random quantum spin model. In 
particular, the influence of randomness on the quantum 
phase transition inherent in the nonrandom spin model 
can be studied. In the presence of randomness the quan- 
tum phase transition becomes a crossover. Our solution 
presented below is based on the Jordan- Wigner transfor- 
mation of the spin Hamiltonian to the Hamiltonian of 
a tight-binding chain of spinless fermions with nearest- 
neighbor and next-nearest-neighbor hoppings and ran- 
dom (Lorentzian) on-site energy. Next-nearest-neighbor 
hopping is a new feature emerging owing to three-site 
interactions which makes further calculations more in- 
volved. We introduce Green functions and flnd exactly 
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the random-averaged Green functions which yield the 
random-averaged density of states. We use the obtained 
density of states to discuss some ground-state and finite- 
temperature properties of the spin model. Although the 
random-averaged density of states can be obtained only 
for a specific probability distribution, the moments of the 
density of states can be obtained for an arbitrary inho- 
mogencous spin chain. These quantities can illustrate 
some general effects on the properties of the quantum 
spin chain caused by inhomogeneity and yields thermo- 
dynamic quantities in the high-temperature limit. Our 
exact results allow to illustrate effect of randomness on a 
quantum phase transition. In particular we discuss how 
a quantum critical region may be modified owing to ran- 
domness. 

The paper is organized as follows. In the next sections 
we define the spin model under consideration (Sec. |TT| 
and calculate the random-averaged density of states 
which yields thermodynamic quantities (Sec. IIII) ). Then, 
in Sec. |lVl we discuss some properties of the spin model 
at zero and nonzero temperatures. We illustrate the ef- 
fect of the introduced disorder for the transverse magneti- 
zation, specific heat, and static transverse susceptibility, 
and put our discussion in a general context of a theory of 
quantum phase transitions with randomness. In Sec. |V] 
we discuss some global properties of the density of states 
for an arbitrary inhomogeneous spin- 1/2 transverse XX 
chain with three-site interactions. Finally, in Sec. IVI[ we 
summarize our findings. 



II. THE MODEL 

To be specific, we consider a linear chain of N spins 
with spin quantum number s = 1/2. Each spin inter- 
acts with spins on nearest-neighboring sites and on next- 
nearest-neighboring sites. Moreover, all spins interact 
with an external magnetic field which acquires a random 
value on each site. The Hamiltonian of the model reads 



-K\ 



H 



«+i4+2)] +E^"<' (2-1) 



where periodic boundary conditions are implied for con- 
venience. Here J and K are the two-site isotropic XY 
(i.e., XX) interaction and the three-site XZX + YZY 
interaction, respectively, and n„ is the transverse field on 
the site n. Although, exact solvability is the main mo- 
tivation to consider the three-site interactions, we note 
that Hamiltonians similar to Eq. (|2.ip may be generated 
in optical lattices<^ 

The on-site transverse fields are assumed to be inde- 
pendent random variables each with the Lorentzian prob- 
ability distribution 
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^ (r!„ - r!o)2 + r2 ' 



(2.2) 



where fig is the mean value and F controls the strength of 
disorder. We are interested in (random-averaged) ther- 
modynamic quantities of the spin model (|2.ip . (|2.2|) . 

As the first step in the calculation of thermodynamic 
quantities of the spin model we perform the Jordan- 
Wigncr fermionization^ to transform the Hamiltonian 
(|2.ip into a bilinear Fermi form 



K f ^ t 



N N ^ N 

~ ^ ] ^ ] CnAnmCm ~ ^ ^ ] ^n- (2-3) 



As typical for the fcrmionic representation of the spin 
model the magnetic field (here its uniform part flo) plays 
the role of a chemical potential. From Ref. i we know 
that the bilinear form in Eq. (j2.3p can be diagonalized. 
After performing the linear canonical transformation 



JV 



N 



n=l 
N 



we find 



1=1 

N N 
i— 1 M— 1 



^ / 1 

^ = E^- r'i'^-" 9 



(2.4) 



(2.5) 



Although this can be done in principle, to find g^n and 
Ay is a complicated task in practice because of nonho- 
mogeneous values of r2„. 

Before we present the solution of the random model, 
for convenience we illustrate briefiy the basic features 
of the nonrandom model, i.e., Qj. = fio is independent 
of the site index n, see Refs. [I6II25I and Fig. |T] The 
nonrandom model exhibits three phases. In the K-VIq 
plane (| J| = 1), the spin-liquid I phase occurs in the re- 
gion — 1 + K/2 < Hq < I + K/2 (dark gray region in 
Fig. [T|), the spin- liquid II phase occurs in the regions 
K < -1/2, 1 + K/2 < Oo < -K/2 - l/i^K) and 
K > 1/2, ~K/2 - 1/(4A') < < -1 + A72 (light 
gray regions in Fig. [1]), whereas the remaining regions 
(white in Fig. [1} correspond to the ferromagnetic phase. 
The two different spin-liquid phases correspond to gap- 
less spinless-fermion systems having two or four Fermi 
points, whereas the ferromagnetic phase corresponds to 
a gapped spinless-fermion system. Although the choice 
of the order parameter which in a transparent way would 
be associated with the modification of the Fermi-surface 
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FIG. 1; The ground-state phase diagram of nonrandom 
model (|2.1 [) with J = ±1 and fl„ = f^o discussed earlier in 
Refs. [I^[2^. Dark gray region corresponds to the spin-liquid 
I phase (two Fermi points), light gray regions correspond to 
the spin- liquid II phase (four Fermi points), and white regions 
correspond to the ferromagnetic phase. 



topology is still under debate r^'i^'^'^ there is no doubt 
that the different phases and the transitions between 
them may be identified looking at the behavior of the 
ground-state transverse magnetization niz as a function 
of Hq or K. Several cusps in the magnetization curve 
indicate quantum phase transition points, see curves for 
r = in Figs. [3] and H below. 



III. THE AVERAGED DENSITY OF STATES 
AND THERMODYNAMIC QUANTITIES 

The frec-fcrmion representation of the spin model, Eq. 
(j2.5p , immediately implies simple formulas for thermody- 
namic quantities, as 



2TJ ' 



f = -T j dLop{Lo) In {2 cosh ; 

P{^) = ^Y.^{^ - (3.1) 



for the Helmholtz free energy (per site). Here we have 
introduced the density of states p{io). The random- 
averaged Helmholtz free energy / is given by Eq. p.ip 



with the random-averaged density of states p{oj), where 



J 



(...) — nn=i / ^^nP{^n){- ■ ■) ■ Thus our task is to find 
p{uj). 

Using p.ip , p.Sp , ()2.4p one can easily convince oneself 
that 



p{uj) 



where 



1 ^ 

T]^E»Gj/..±ie), 



(3.2) 



1 

2^ 



G',T,„(t) = Ti0(±t)({c„(t),4}), 

dw exp(— iti;t)Gj„^(aj ± ie) (3.3) 



[9{x) is the Heaviside step function] are the retarded and 
advanced temperature double-time Green functionsi^"— 
On the other hand, one easily finds the following set of 
equations for GJf,„ (w ± ie) 



J 



{uj±ie-nn)Gf^^(uj±\e) 
[G^i,™(^±ie) + Gj+i,„(a.±ie)] 



+ ^ [G,T-2,m(^ ± ie) + G,T+2,,n(^ ± i^)] = '5„™. (3.4) 

Because of nonhomogeneous values of 51„ it is not possi- 
ble to solve p.4p and to find the required diagonal Green 
functions G^„(a; ± ie) which enter Eq. p.2[) . However, 
it is well knowriF^ that if f2„ is a Lorentzian random 
variable (|2.2p the set of equations p.4p can be averaged 
over random realizations leading to a set of equation for 
translational-invariant random-averaged Green functions 

Gnm(aj). Supposing that ^ln is a complex variable and 
noticing that G,7,„(w+ie) [G^,„(w— ie)] cannot have a pole 
in the lower [upper] half-plane of the complex variable 
we perform the averaging with (j2.2p by means of contour 
integrals closing the contours of integrations in the half- 
planes where the Green function has no polesJ^"— As a 
result we obtain 



(w ± ir - r^o) Gtm{oj) 



K 



G',T-l,™(^) + Gj+l,m('^) 



(3.5) 



The set of equations p.Sp possesses translational symmetry already and therefore 

exp [i(n — m)K\ 



1 

G5fm(w) = TT I dK- 

J-TT 



Lo ~ VLq ~ J cos K + -|- cos(2k) ± ir 
To evaluate the integral in p.6p we introduce a new variable z ~ exp(iK). Then Eq. (|3.6p becomes 

27ri 



Gjm(w) 



yii — m+l 



(z4 + 1) - lz{z^ + 1) + (w - r!o ± ir)z2 ' 



(3.6) 



(3.7) 
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where the contour of integration runs counterclockwise alon g the unit circle in the complex plane z. 



The calculation of (|3.7p is simple if cither X = 
or J = yielding G,T„(w) = l/^(w - f^o ± ir)^ - J2 

(see Ref. [13) or G^F„(a;) = l/^(a; - f^o ± ir)^ - X2/4. 
For arbitrary values of the interaction constants, < 
\KlJ\ < oo, we have to solve the 4th order algebraic 
equation 



K K K 



(3.8) 



which is a quasi-symmetric one (i.e., of the form aoz* + 
aiz^ + 02^^ + aimz + aom? = with m = 1). Dividing 
Eq. p.8p by and using the variable change y — z + 1/ z 
we immediately find 

J 



2/± = ±5, 



if 



(w - r^o ± ir) + 2. 



As a result, 



1, 



(3.9) 



(3.10) 



where y is either y+ or Let us denote the roots of 
Eq. (1331), which are given in Eqs. ((XTU)) and ((^ . by 
zi, 22, 2:3, Z4, |zi| < \z2\ < \z3\ < |z4|, scc Appendix [J] 
Only two roots are inside the unit circle l^^l < 1 resulting 
in 
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K 



Zl 



[zi - Z2) {zi - 23) (zi - Z4) 

Z2 

{Z2 - Zl) (Z2 - Z3) (Z2 - Zi) 



(3.11) 



Then the density of states piu)) is calculated according 
to Eq. p.2p . The described scheme for F = reproduces 
the density of states of the nonrandom model reported 

in Ref. [H 

Our results for p(a;) for typical sets of parameters are 
shown in Fig. [21 We put VLq ~ 0, since a nonzero fio 
leads only to a trivial shift along the w-axis, see, e.g., 
Eq. p.6p . At the van Hove singularities present for the 
nonrandom model in Fig.|2lthe density of states exhibits 
the typical one-dimensional inverse square-root singular- 
ity. For parameters fio or K where a quantum phase 
transition occurs, a van Hove singularity is located at 
w = fio (as in Fig. |2}d, where it is at cj = fio = 0). 
The middle peak in the density of states shown in Fig. [21 
(it appears if \K\ > 1/2, see Fig. [1} indicates a quantum 
phase transition between the spin-liquid I and spin-liquid 
II phases. Small randomness leads mainly to rounding of 
the van Hove singularities and to the appearance of tails 
in the density of states above and below the band edges 
of the nonrandom model. With increasing of F the den- 
sity of states becomes more and more smeared out, i.e.. 
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FIG. 2: The random-averaged density of states p{uj) for the 
spin model ((23), with J ^ 1, K = 0.75 (a), K ^ 2 

(b), Qo = 0, and r = (solid), T = 0.05 (dashed), P = 0.1 
(dotted). 



the fingerprints of the van Hove singularities disappear 
and the tails increase. These changes of the density of 
states due to randomness influence the ground-state and 
finite-temperature properties to be discussed in the next 
section. 



Finally, it is important to note the following. A ground 
state of the nonrandom spin model (|2.ip is generally 
speaking a product of one-particle (more precisely, one- 
fermion) states, which arc plane waves. While for F = 
the one-particle states within the band(s) are extended 
states with infinite localization length, for any small (di- 
agonal) Lorentzian disorder, i.e., for F > 0, all one- 
particle states of the Hamiltonian (|2.5p become localized, 
see, e.g., Ref. Isol As a result, the nature of the quantum 
phases of the nonrandom model changes and the quan- 
tum phase transition inherent in the nonrandom model 
becomes a crossover 1^ 
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IV. GROUND-STATE AND THERMODYNAMIC 
PROPERTIES 

The obtained (random-averaged) density of states per- 
mits to examine various quantities characterizing behav- 
ior of the spin model at zero and nonzero temperatures, 
see Eq. ()3.ip . For the ground-state energy, the entropy, 
and the specific heat we have 



(a) 



= / du!p{u 



Co = y"da;p(a;)Jy, 

in 2 cosh — — tanh — - 

V 2TJ 2T 2T 



I 2T 

\cosh 5y 



(4.1) 



(4.2) 



(4.3) 



respectively. Next, for the transverse magnetization and 
the static transverse susceptibility we have 



_ df 
olio 



dm. 



1 

"2 

1 

"4T 



du!p{uj) tanh 



2r' 



dujp{uj) 



cosh ^ 



(4.4) 



(4.5) 



respectively. 

In Figs. [3] and |4] we report some results for the ground- 
state transverse magnetization. In the ground state of 
the nonrandom spin system the model can be in three 
different phases (spin-liquid I, II, and ferromagnetic), see 
Rcfs. llqKSi and the discussion in Sec. HT] The transitions 
between them can clearly be detected by the cusps in the 
magnetization curves in Figs. [3]and|4l 

Let us briefly discuss a prominent feature of the mag- 
netization curve, namely the steep part in the curve near 
saturation {Hq ~ -0.708) seen in Fig.[5]3 for F = 0. This 
jump-like behavior resembles the magnetization jumps 
observed in frustrated quantum antiferromagnetsj^ The 
corresponding density of states, sec Fig. [5^, shows a nar- 
row upper band, present for K > 1/2. The two singu- 
larities defining the band edges approach each other if 
K — > 1/2, i.e., the upper band becomes flat. However, 
by contrast to flat bands discussed in Ref. [s^ the number 
of states in the narrow upper band of our model decreases 
with decreasing of band width. As a result the middle 
cusp in the magnetization curve, see Fig. [3]d, related to 
the middle singularity in the density of states, see Fig. [5^, 
moves to the left cusp this way yielding the steep part be- 
fore saturation seen in Fig. [3Jd. The slope of that part of 
the magnetization curve increases if K ^ 1/2, however, 
at the same time its height decreases and vanishes finally 
at if = 1/2, where the magnetization approaches the sat- 
uration continuously with an infinite slope, see Fig. 
Moreover, 1/2-m^ cx (Jig -1-0.75)^ if f^o+0. 75 4-0 with 
e = 1/4 instead of the usual value e = 1/2, see Ref. [2]. 
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FIG. 3: Ground-state transverse magnetization — (|4.4p 
versus Q,o for the spin model 1)2. II) . (|2.2p with J = 1, K = 
0.5 (a), K = 0.75 (b), K = 2 (c) for T = 0, ...,0.2. The 
quantum phase transitions occurring in the nonrandom model 
at Qo = -0.75, fio = 1.25 (a), at Qo « -0.708, fio = -0.625, 
Q.0 = 1.375 (b), and at ilo = -1.125, f2o = 0, f2o = 2 (c) are 
signaled by kinks in the magnetization. 



The effect of randomness on the magnetization mz is 
similar to that of a finite temperature. For small random- 
ness the cusps in the mj-curvcs, which indicate bound- 
aries of different ground-state phases, become rounded, 
indicating that a quantum phase transition present at 
r = transforms into a crossover at F > 0. Although 
even small randomness is sufficient to erode the bound- 
aries between different ground-state phases by a notice- 
ably rounding of the cusps of mj, it may have almost 
no influence on for parameter values corresponding 
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K 

FIG. 4: Ground-state transverse magnetization —ml (|4.4p 
versus K for the spin model (|2.1I) . (|2.2I) with J = 1, fio = 
(a), fio = 0.5 (b), ^0 = 0.75 (c) for T = 0, ...,0.2. The 
quantum phase transitions occurring in the nonrandom model 
at K = -2, K = 2 (a), at K = -1, K = 3 (b), and at 
K = —1, K — —0.5, K = 3.5 (c) are signaled by kinks in the 
magnetization. 

to the spin- liquid phases, see Figs. [3] and ID Other pecu- 
liarities of the nonrandom model, namely nonzero mag- 
netization at r^o = 0; zero magnetization at nonzero fig, 
as well as saturated magnetization for |f7o| < \J\; see 
Fig. [31 become less pronounced as the strength of disor- 
der increases. 

In Figs. [5] and [6] we report some of our findings for 
nonzero temperatures (for the sake of brevity we consider 
the case fio = 0, only). The specific heat c for various 
parameter sets is presented in Fig. [51 where we show c/T 
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FIG. 5: c/T versus T, c is the specific heat given in Eq. (|4.3|l . 
for the spin model ^J}, with J = 1, A' = 1, 1.9, 2, 2.1 
(from top to bottom), Slo = 0, and F = 0, . . . , 0.2. 
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FIG. 6: Static transverse susceptibility — Xzz (|4.5p versus T 
for the spin model ^J}, ([Ijj with J = 1, A' = 1, 1.9, 2, 2.1 
(from top to bottom), Slo ~ 0, and F = 0, . . . , 0.2. 



as a function of temperature T. For the nonrandom case 
F = it is knowni^ that c(T) oc T in the limit T ^ for 
all K except K = i^Tcrit = i2|J|. For the critical value 
of K, K ^ KcTit, we have c(T) oc vT. 

Let us start with the discussion of the behavior oic/T 
for K = 1, i.e., the system is quite far away from the 
quantum critical point K^it = 2. Then we have the typ- 
ical high-temperature maximum in c{T) related to the 
relevant energy scale (depending on J and K) and a con- 
stant value of c/T at low T corresponding to c cx T. 
The position of the high-temperature maximum moves to 
higher temperatures as K increases. The overall modifi- 
cation of the c/T versus T curve by randomness is small. 
There is only a small change of the slope in the depen- 
dence c on T as T — )■ and a slight shift of the height 
and the position of the high-temperature maximum. At 
the quantum critical point K = i^crit = 2 we have a com- 
pletely different behavior of c/T. The high-temperature 
maximum is still present (and again there is only a weak 
effect of randomness on that maximum). However, below 
the maximum there is an increase of c/T with decreasing 
T indicating the c oc Vt dependence. While this in- 
crease is monotonous till T — > for the nonrandom case, 
in the random system there is only a finite region of T 
where this increase of c/T can be observed. At lower T 
again the c oc T regime sets in, i.e., the randomness de- 
stroys the -s/T-dependence in favor of the T-dependence 
as T ^ 0. A similar behavior can be found for K values 
near i^crit = 2, e.g., at K = 1.9 or 2.1, i.e., we observe 
the A/T-dependence for the low-temperature specific heat 
in a finite temperature range below the high-temperature 
maximum. That is the typical quantum critical behavior 
appearing in the vicinity of a quantum critical pointii"— 
Interestingly, in the specific model under consideration 
for small randomness there is a second maximum for c/T 
(but not for c) at a lower temperature T* (which is a 
reminiscence of the singularity for K ~ A'crit at T = 0) 
before the c (xT regime sets in at very low temperatures. 
According to the above discussion for Figs. Eh, Et, and 
EJi, we conclude that for the system at = i^crit at low 
but finite temperatures the randomness has a similar ef- 
fect as shifting the nonrandom system slightly away from 
the quantum critical point. Note that further increasing 
of F removes the low-temperature maximum in c/T and 
also the critical-like VT behavior of c disappears, see the 
curves for F = 0.1 in panels for K = 1.9, 2.1 of Fig. [5l 

The temperature dependence of the static transverse 
susceptibility Yzz shown in Fig. [6] exhibits many similar- 
ities to that of c/T discussed above. In the nonrandom 
case it is knowrnl^ that for K = Kent one has Xzz oc 1/VT 
(critical behavior) as T — > 0, whereas Xzz remains finite 
at T = for noncritical values of K. For K around /STcrit 
a reminiscent of the critical behavior emerges in the low- 
temperature region starting from certain finite tempera- 
tures. For small nonzero F, Xzz may exhibit the critical 
behavior in a certain temperature range starting from fi- 
nite temperatures if K is equal to or is close to i^critj 
see the curves for F = 0.01 in panels K = 1.9, 2, 2.1 of 
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Fig.l 

We can use the above discussed observations of the 
temperature profiles of c/T and Xzz to construct a phase 
diagram of the random quantum spin chain (j2.ip , (j2.2l) , 
i.e., to determine the quantum critical region in the K-T 
half-plane, see Fig. [71 As an indicator of the critical (i.e., 
\/r-likc) behavior we choose the value of the derivative 
—d{c/T)/dT. To be more specific, we use the circum- 
stance that for K in the vicinity of i^crit and small F > 
c/T exhibits a maximum at T* (for K = Kcmt, F = it 
is a divergency at T = 0) and a l/\/T-\\ke decrease (due 
to c cx a/T) in a certain temperature region above T*. 
As discussed above this behavior is a trace of the quan- 
tum critical point. The plots shown in Fig. [7] are based 
on a quantitative analysis of —d{c/T)/dT and show the 
value of —d{c/T)/dT as grayscale plots. All white areas 
in this figure belong to negative values of —d{c/T)/dT. 
The lower boundary of the quantum critical region is re- 
lated to T* (the temperature where the low-temperature 
maximum of c/T is located), since for T < T* one has 
—d{c/T)/dT < 0, whereas already for T slightly above 
T* the derivative —d{c/T)/dT > becomes quite large. 
As a result, the lower boundary of the quantum criti- 
cal region is quite sharp in Fig. [71 Note that for F — 
our plot reproduces the typical picture for a quantum 
critical region.—"— For small F > the lower boundaries 
in the K-T half-plane remain nearly straight lines with 
different slopes « —0.225 and sa 0.096 below and above 
K = 2, respectively. We compare these lower bound- 
aries for various strengths of disorder F in Fig. [8l where 
the numerically determined values of T* as a function 
of K are shown. It is obvious that a noticeable change 
in the slope is observed only in a small vicinity of the 
K = 2. On the other hand, the upper boundary of the 
quantum critical region is not sharp, since ~d{c/T)/dT 
varies smoothly, if the temperature is further growing. 
From Fig. [7] it is obvious that for small randomness (see 
the panel for F = 0.01 in Fig. [3 and the corresponding 
curves in Fig. [5]) at low temperatures the area, where sig- 
natures of quantum critical behavior can be observed, is 
enlarged, whereas a further increase of randomness then 
leads to a shrinking of that area. For F > 0.1 the signa- 
tures of quantum criticality completely disappear. 

We have to say in the end, that the adopted criterion 
to indicate the effect of randomness on quantum critical 
behavior represents only one possibility. In fact, with 
this criterion we only indicate a region where the specific 
heat c behaves like T" with a < 1 and claim that such a 
behavior is a remnant of the quantum critical one in the 
nonrandom model when a = 1/2. 
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V. GENERALIZATION OF THE MODEL 

Our consideration until now was restricted to the XX 
two-site interaction and the XZX + YZY three-site in- 
teraction [see Eq. p.ip ] with a Lorentzian probability 
distribution of the random transverse field, which allows 



FIG. 7: Critical region as indicated by the value of 
~d{c/T)/dT > in the half-plane K-T for the spin model 
dSU, with J = 1, no = and r = 0, 0.01, 0.05, 0.1 
(from top to bottom). See further explanations in the main 
text. 
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FIG. 8: The lower boundary of the quantum critical region 
determined by the position of the maximum T* of c/T for 
the spin model (pTTj) . (f2^ with J = 1, fio = and T = 
0, 0.01, 0.05. 



the exact calculation of averaged density of states. In this 
section we want to generalize the model by (i) extending 
the Hamiltonian and (ii) considering general probability 
distributions. 

First we illustrate the extension of the Hamiltonian. 
We can add to the Hamiltonian (|2.ip a Dzyaloshinskii- 
Moriya two-site interaction with the strength D„ and a 
X ZY — Y Z X three-site interaction with the strength i?„, 

n 

<r)(„x„v — 1 

+En (SnSra+lSra+2 ~ + ^ ^nSn- (5-1) 

n 

In fcrmionic representation Eq. (j5.ip reads 

H An ( 4c« ^ ^ ) + {Bncicn+l + H.C.) 



-f (C„c];c„+2+H.c.)] , (5.2) 

where An = ^n, Bn = (J„ -I- i£'„)/2, and C„ = -(A'„ + 
ii?„)/4. Comparing Eq. (|5.2p with Eq. (|2.3p we con- 
clude that for a Lorentzian probability of r2„ and uni- 
form parameters Jn = J, Dn = D, Kn = K, and 
En = E the formula (|3.6p for the averaged Green func- 
tion is valid also for the generalized model (|5.ip if we 
change Jcosk — (AY2)cos(2k) — >■ Jcosk + DsinK — 
{K/2) cos{2k) — {E /2) sm{2K). Further analysis can be 
performed as in the previous sections. Although the re- 
sults obtained from the exact averaged density of states 
are more general now, they do not give basically new 
features in comparison to those discussed above. 

So far our discussion has been referred to a special 
kind of randomness, i.e., to a Lorentzian transverse field. 
Only in this case we can provide an exact analysis of 
the thermodynamics of the spin model. However, it is 



possible to analyze some global properties of the averaged 
density of states (which imply corresponding properties 
of the spin chain) for an arbitrary inhomogeneous spin- 
1/2 XX chain with three-site interactions in a transverse 
field (|5.ip . For that we consider the moments of the 
density of states p.l|> ^ 



n=l 
1 ^ 

du:Lop{^) = -J2{{[cn,H],ci}), 

n=l 
1 ^ 

dcouj^picj) = - Y,{{[[cn,H],H],ci}) , 

71=1 
1 ^ 

dLou^piu) = - ^ ({ [[[c„, H],H],H], ct }) 



(5.3) 



etc. Calculating the right-hand sides in Eq. (|5.3p with 
the Hamiltonian (15.21) we find 



n 



n 

+Al + 2 {\Cn-2 
|2 



\Bn-l\^ + \B„ 



\Bn-l\^An-l 



\Cn 



\Bn-l\'' + \Bn\^ + \Cn\^)An 



\Bn\ An+l + \Cn\ An+2 
25R {C:^^Bn-2Bn-l + C:_,Bn-lBn + QBnBn+l)] 



(5.4) 



The moments of the random-averaged density of states 
p{uj) follow from Eq. (|5.4p after a corresponding averag- 
ing. 

We may use the derived moments of the density of 
states (|5.4p to examine some general properties of the 
(homogeneous or inhomogeneous) spin chain (j5.ip . For 
example, the ground-state transverse magnetization is 
given by the formula 



dwp(cj) 



h4 



(5.5) 



see Eq. (|4.4p . It has been shown that the uniform spin 
model (j5.ip may exhibit a nonzero transverse magneti- 
zation niz in zero transverse field r2„ = fig = 0'^ ^^^''^^ 
From Eq. (|5.5p it is clear that uiz 7^ at zero field if the 
density of states is asymmetric, i.e., if the third moment 
of the density of states at zero field is nonzero. 

In the uniform (nonrandom) case we have M(3) = 
A^ + 6(|e|2 + \C\^)A + 6'iSi{C*B^i or in the zero-field case 
M(3) = 6n{C*B^) = -{3/8)[{J^ - D'^)K + 2JDE]. From 
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the latter formula it is obvious that only three-site inter- 
actions may lead to nonzero magnetization. More pre- 
cisely, the XZX + YZY interaction if ^ D"^ or the 
XZY - YZX interaction if JL* 7^ leads to 7^ at 
zero field. Formulas for the moments of the density of 
states (|5.4|) permit to examine the effect of randomness 
on this property. An example of such analysis is given 
in Ref. Il5b where some consequences of the correlated 
off-diagonal and diagonal disorder were discussed. 

Finally we note that the high-temperature properties 
of the spin model are determined by the lower moments 
of the density of states [see Eq. p.ip ] and therefore the 
thermodynamic quantities in the high-temperature limit 
can be examined accurately for any type of disorder on 
the basis of Eq. (fO]) . 

VI. CONCLUSIONS 

We have considered a random spin- 1/2 XX chain with 
three-site interactions. For random on-site transverse 
field with Lorentzian probability distribution we have 
calculated exactly the random-averaged density of states 
and the corresponding random-averaged thermodynamic 
quantities of the model. For arbitrary inhomogeneous 
Hamiltonian parameters we have calculated the three 
first moments of the density of states which determine 
some general properties of the spin model and yield its 
thermodynamic quantities in the high-temperature limit. 
The effect of a random transverse field on the ground- 
state magnetization process and on the temperature be- 
havior of the specific heat and static transverse suscep- 
tibility has been analyzed. As a main result we have 
discussed how the quantum critical behavior of the spin- 
1/2 XX chain with three-site interactions is modified by 
randomness. While for large enough randomness all sig- 
natures of quantum critical behavior disappear, we find 
even a slightly enlarged temperature area for very small 
randomness, where such signatures can be observed. 

Finally, we have to underline that our analytical find- 
ings are restricted to the random-averaged density of 
states and thermodynamic quantities. To our best knowl- 
edge, almost all applications of Lloyd's model are re- 
stricted to one-particle quantities which can be obtained 



rigorously. The calculation of the two-particle quantities 
meets notorious difficulties and requires approximations, 
see, e.g., Ref. [ssl On the other hand, other quantities of 
interest which are not related to the one-particle Green 
functions (j3.3p but are important for following final out- 
comes of introduced randomness (e.g., two-spin corre- 
lation functions) may be examined only numerically for 
finite chains?^ see also Ref.H^- However, numerical stud- 
ies are beyond the scope of the present paper and will be 
reported separately. 
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Appendix A: The roots of Eq. (|3.8|l 

In this appendix we give explicit expressions for the 
roots of Eq. (123). They read: 



{2:1,^2, 2:3, -24} = 



K 


-.9- ^ 




-9, 
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- 4 
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+ ,g - 1 


/(^ 


+.9; 
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- 4 
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- ,g + 1 




-a, 
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- 4 
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f ,9+ ■\ 


hi- 
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-4 



(Al) 



with g = ^{J/Kf - {^/K){uj - r^o + iF) + 2 in the 
case of the retarded Green functions p.3p or g = 
y/{J/Kf - {'ilK){uj - r^o - ir) + 2 in the case of the 
advanced Green functions (|3.3p . 
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